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We investigate a formulation of Poincare invariant quantum mechanics where the dynamical input 
is Euclidean invariant Green functions or their generating functional. We argue that within this 
framework it is possible to calculate scattering observables, binding energies, and perform finite 
Poincare transformations without using any analytic continuation. We demonstrate, using a toy 
model, how matrix elements of e~^ H in normalizable states can be used to compute transition 
matrix elements for energies up to 2 GeV. We discuss some open problems. 

o": 

We investigate the possibility of formulating Poincare invariant quantum models of few-body systems where the 
dynamical input is given by a set of Euclidean-invariant Green functions. This is an alternative to the direct construc- 
. tion of Poincare Lie algebras on few-body Hilbert spaces. In the proposed framework all calculations are performed 
using time Euclidean variables, with no analytic continuation. 

One potential advantage of the Euclidean approach is that it has a more direct relation to Lagrangian based field 
theory models. One of the challenges is the construction of a robust class of suitable model Green functions. In this 
paper do not address this problem; wc assume that this has already been solved and discuss how one can calculate 
£C) • observables without analytic continuation. 

Most of what we propose in not new, it is motivated the reconstruction theorem of a quantum theory in Euclidean 
field theory. The fundamental work was done by Osterwalder and Schrader [l[ Q . The approach illustrated in this 
work is strongly motivated by Frohlich's Q elegant solution of the reconstruction problem using generating functionals. 
One of the interesting observations of Osterwalder and Schrader is that locality is not needed to construct the quantum 
' theory. 

To keep our discussion as simple as possible wc assume that we are given a Euclidean invariant generating functional 
for a scalar field theory. This input replaces the model Hamiltonain. We assume that this generating functional is 
Euclidean invariant, positive, reflection positive, and satisfies space-like cluster properties. These requirements are 
defined below. The conditions on the generating functional imply conditions on Green functions in models based on 
a subsets of Green functions. 

For a scalar field the Euclidean generating functional Z[f] is the functional Fourier transform of the Euclidean path 
measure: 



< 



43 
Oh 



> 
00 

O 

" f D Ule-^M+W) ^ (i) n 



oo 

o 

O 



fD e [<f>]e- 

n times 



where /(x) = /(r, x) is a test function in four Euclidean space-time variables and S n (xi,--- ,x„) is the n-point 
Euclidean Green function. 

The generating functional is Euclidean invariant if Z[f] = Z[f] where /'(x) = f(E^ 1 (x — a)) where x — > Ex + a is 
Jv>( \ a four-dimensional Euclidean transformation of the arguments of /. 

The generating functional is positive if for every finite sequence of real test functions {/j} the matrices Eij = 
' Z[fi — fj] are non-negative. 

The generating functional is reflection positive if for every sequence {fi}, of real test functions with support for 
positive Euclidean time, the matrices My = Z[fi — Qfj] are non-negative, where (0/)(r, x) = /(— r, x) is Euclidean 
time reflection. 

The generating functional satisfies space-like cluster properties if 

lira (Z[f + g a ] - Z[f]Z[g}) ^ (0.2) 

|a| — >oo 

where 

.9a (T, x) := g(r, x - a). (0.3) 



These are the primary requirements that are expected of an acceptable generating functional. 
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I. HILBERT SPACE 



We begin by representing vectors by wave Junctionals of the form 

N b N a 

B[4>] = bjeWA C[<t>] = J2 c k ^ 4(9k) (1.1) 

3=1 k=l 

where bj and Cfc are complex constants and /j(x) and <?fc(x) are real Euclidean test functions. The argument "0" plays 
the role of a formal integration variable. 

We define a Euclidean-invariant scalar product of two-wave functionals by 

N b N c 
3 = 1 k=l 

This becomes a Hilbert space inner product by identifying vectors whose difference has zero norm and adding con- 
vergent sequences of finite sums. We call this space the Euclidean Hilbert space. 

Reflection positivity can be used to define a second Hilbert space. Vectors arc represented by wave functionals of 
the form (|1 . ip where the test functions /j(x), gfc(x) are restricted to have support for positive Euclidean times. We call 
these test functions positive-time test functions. We define the physical scalar product of two such wave functionals 
by 

N b N c 
j=l k=l 

As in the Euclidean case, this becomes a Hilbert space inner product by identifying vectors whose difference has zero 
norm and adding convergent sequences of finite sums. We will refer to the resulting Hilbert space as the physical 
Hilbert space. Reflection positivity is equivalent to the requirement that 

(B\B) > 0. (1.4) 



II. POINCARE LIE ALGEBRA 

Note that the determinant of the 2x2 matrices 

x= ft-zx-iy\ x= fir-zx-iy\ ^ _ 

\x + iyt + zj I x + iy it + z J \ ■ ) 

gives the Lorentz and Euclidean invariant distances. The determinants are preserved under the linear transformations 
X — > X' = AXB 1 and X — > X' = AXB* where A and B are complex matrices with determinant 1. These transfor- 
mations are generally complex but the determinants remain real. It follows that the pair (A, B) equivalently defines 
both complex Lorentz and complex 0(4) transformations. Real Lorentz transformations have B = A* while real 
0(4) transformations have A and B € SU(2). In this section we use the observation that real 0(4) transformations 
correspond to complex Lorentz transformations to extract Poincare generators on the physical Hilbert space. 
Finite Euclidean transformations, T(E,a), act on wave functionals as follows 

N„ 

T(E,z)B[ct>\='£b j e i *^\ (2.2) 
j=i 

where /B, a j(x) = fj(E^ 1 (x — a)) and E G 0(4). Since real Euclidean transformations preserve the Euclidean scalar 
product (•, •), T(E 1 a) is unitary on the Euclidean Hilbert space. 

These same transformations, with restrictions on the domains and group parameters to ensure the positive time 
support condition is preserved, are defined on the physical Hilbert space, but the resulting transformations are not 
unitary. 

For three-dimensional Euclidean transformations, T(E,a) maps the physical Hilbert space to the physical Hilbert 
space in a manner that preserves the physical Hilbert space scalar product. This implies that for space translations 
and ordinary rotations T(E, a) is unitary on the physical Hilbert space. 
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Positive Euclidean time translations, T(I,((3,0)), j3 > 0, map the physical Hilbert space to the physical Hilbcrt 
space, however because of the Euclidean time reversal operator in the physical scalar product, Euclidean time transla- 
tions are Hermitian, rather than unitary. It is possible to use the unitarity of on the Euclidean Hilbert space along 
with reflection positivityQ to show that positive Euclidean time evolution is a contractive Hermetian semigroup on 
the physical Hilbert space. 

Rotations in planes that contain the Euclidean time direction do not generally preserve the positive Euclidean time 
support constraint. However, if the test functions are restricted to have support in a cone with axis of symmetry along 
the Euclidean time axis that makes an angle less than n/2 with the time axis, then rotations in space-time planes 
through angles p small enough to leave the cone in the positive-time half plane are defined on this restricted set of wave 
functionals. On this domain and for this restricted set of angles Euclidean space-time rotations are Hermitian. They 
form a local symmetric semigroup. What is relevant is that just like one-parameter unitary groups and contractive 
Hermitian semigroups, local symmetric semigroups have self-adjoint generators Q @0- 

The result is that on the physical Hilbert space the various one-parameter subgroups of the real Euclidean trans- 
formations have the form 

T(E, a) -> e la P , e l6> J , e~ pH , e^ K (2.3) 

where H, P, J, K are all self-adjoint operators on the physical Hilbert space. It can also be shown by direct calculation 
that the infinitesimal generators satisfy the Poincare commutation relations. This is a consequence of the relation 
between the complex Lorentz and complex 0(4) groups. 

Matrix elements of the generators can be computed by differentiating T(E, a) with respect to the group parameters: 

(B\3\C) 



(B\P\C) 



(B\H\C) 



(B\K\C) 

where p is the axis and angle of a rotation in a plane containing the Euclidean time direction (it is an imaginary 
rapidity). 

In this section we have illustrated how the Poincare generators can be constructed directly from the Euclidean 
generating functional without using any analytic continuation. 



g N b N <= 

l d8 £ E fykZlSk - ®f E( 0),oJ\0-o 
j=i fe=i 



(2.4) 



g N b N c 

"ifa £ £ b*c k Z[g k - e/ /iaj ]| a=0 



(2.5) 



3 = 1 k=l 



d(3 



(2.6) 



3 = 1 k=l 



df> j=l k=l 



31 1/3 = 



(2.7) 



III. PARTICLES 



Particles are associated with eigenstates of the mass Casimir operator of the Poincare group with eigenvalues in 
the point spectrum. Matrix elements of the square of the mass operator are 

Q2 Q2^]^^ 



(b\m 2 \c) - ( W2 + p EEW ~ e/w^Wo (3.1) 

J 3 = 1 k=l 



Since the positive-time wave functionals are dense in the physical Hilbert space it is possible to construct an 
orthonormal basis of wave functionals {i?„[0]} satisfying 

B n [<t>] (B n \B m )=S mn . (3.2) 
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Point eigenstates of the mass operator are normalizable solutions of the eigenvalue problem 

((M 2 - X 2 )B X ) [4>] = 0. (3.3) 
In the orthonormal basis {B n [<fi]} this eigenvalue equation becomes 

B x [(f>} =J2 b nB n [4>] ^(B m |Af 2 |B„)fo„-A 2 6 m (3.4) 

n n 

where the sum is generally infinite. 

States of sharp linear momentum and canonical spin can be extracted using translations and rotations. Specifically 
mass-momentum eigenstates are given by 

(C\Bx( P )) = f 7^7je- 4 P- a £^& IlC ;Z[/ / , (0 , a) , l -e. 9 ,] (3.5) 

which can be normalized so 

(B x (p')\B x (p)}=8(p'-p) (3.6) 
Similarly, it is possible find simultaneous eigenstates of mass, linear momentum and spin using 

(C\B Xtj (p,f,)) := [ ^ di?(C|T(i?,0)|i? A (i?- 1 p))^;( J R) (3.7) 

JSU(2) u= _ j 

where dR is the SU (2) Haar measure. 

The wave functional B X j{p, fi)[<f>] describes a particle of mass A, linear momentum p, spin j and z-component of 
canonical spin /i. 

We remark that if this state is non-degenerate, then it must transform irreducibly with respect to the Poincare 
group. This means that 

(C\U(A,a)\B x ^(p,ri) := £ / dp' (C\B Xij (p' , v'))V^,. p JA, a] (3.8) 

I*' =-3 

where 

Vpi'vJ^*} ■= ((A,i),p / J /i , |C7[A J a]|(A,i),p,M) (3.9) 

is the known Wigner fucntion of the Poincare group in the basis |(A, j),p, fj). 

As emphasized in the previous sections, all of the calculations were done using Euclidean Green functions and test 
functions, with no analytic continuation. Equation (|3.8[) demonstrates how to perform finite Poincare transformation 
on the one-body solutions. 



IV. SCATTERING 



The conventional treatment of scattering problems in quantum field theory is formulated using the LSZ asymptotic 
conditions. These have the advantage that they can be implemented without solving the one-body problem, which is 
non-trivial in field theories. However, given one-body solutions it is also possible to formulate scattering asymptotic 
conditions using strong limits. These asymptotic conditions were given by Hagg and Ruclle 8] [9], and they are the 
most natural generalization of the formulation of scattering that is used in non-relativistic quantum mechanics. 

In this work it is useful to use a two Hilbert space formulation of Haag-Ruelle scattering theory [ll| [l2j , where an 
asymptotic Hilbert space is introduced the formulate the asymptotic conditions on the scattering states. All particles 
appear as elementary particles in the asymptotic space; the internal structure (bound-state wave functions) appear 
in the mapping to the physical Hilbert space. 

For a scalar field theory with a mass eigenstate with eigenvalue A, Haag and Ruelle multiply the Fourier transform 
of the field by a smooth function p x (p 2 ) that is 1 when p 2 = —A 2 and vanishes when —p 2 is in rest of the mass 
spectrum of the system. The product, 4> p {p) := (f>(p)p x (p 2 ), is Fourier transformed back to configuration space. The 
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resulting field, <fi p (x), has the property that it creates a one-body state of mass A out of the vacuum. It transforms 
covariantly, but is no longer local. While this is not a free field, it asymptotically looks like a free field, and it is useful 
to extract the linear combination of 4> p (x) and (j)p( x ) that asymptotically becomes the creation part of the field: 

A(f,t) := -i J cf> p (x) d a f(x)dx (4.1) 

where f(x) is a positive-energy solution of the Klein Gordon equation with mass A. Haag and Ruelle prove that the 
scattering states of the theory are given by the limits: 

lim |||*±(/i > ---/n))-A(/ n> t)-..A(/ 1> t)|0)|| = 0. (4.2) 

t— >±oo 

To express this in a two Hilbcrt space notation we rewrite A(f, t) as 

A(/,f)|0) = e im J ([H,4> p (p)]-u, x (p)4p(pj) e-' iffi |0)dpe-^(P)*/(p)dp. (4.3) 

It follows that 

A(f n> t) ■ ■ ■ A(f u t)\0) = e im I ^(p„)---A(p 1 )|0)dp„---d Pl e- iff »74Pn)---/i(Pi) (4-4) 

J S , ' 

where 

ff = £> A ( Pi ). (4.5) 

3 

In this notation equation (|4.20 has the form 

|*±(/i, •••/»)> = I™ e iHt $e- iHat \f) = fi±|f). (4.6) 

t— >oo 

This can be expressed in the Euclidean generating functional representation by replacing <p(p)p\{p 2 ) by the B\j(jp, /i) 
which also creates a one-body state of mass A out of the vacuum. The operator $ becomes 

$(p„,/i„---pi,^i )[</>] := \J[([H,Bx h , jk (pk^k)] ~ u\ k (Pk)B\ k , jk (p k , Mfc))) [<t>] (4-7) 
where the wave functionals are treated as multiplication operators and 

Two Hilbert space wave operators are defined by 

|*±(/i, "•/»)> = I™ e iHt $e- iHot \f) = n±|f>. (4.9) 

t— >oo 

The wave operators satisfy 

l7[A,a]fi± = fi±E7>[A,a] where U f [A,a] = <g>Ux k , jk [A,a}. (4.10) 

Since the asymptotic particles transform like free particles with physical masses, this formula can used to compute 
finite Poincare transforms of scattering states. 

V. COMPUTATIONAL CONSIDERATIONS 

One of the difficulties with using the generating functional representation to do scattering calculations is that we 
have no simple means to construct e lHt on the physical Hilbert space. This can be overcome using a trick. In 
non-relativistic scattering theory Kato and Birman [ll[ [l2j showed that if 

lim e lHt $e-' lHot \ij) = \ij)±) (5.1) 

t— »±oo 
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then for admissible functions x 

lim e <x ^*$e- <x ^ ^|^) = \ip±). (5.2) 

t— >-±oo 

A useful choice of an admissible x is x( x ) — —e~^ x for f3 > 0. If this result is used in (|5.2j) then we have the alternative 
representation of the scattering state 

|V±)= lim e- me ~" H <S>e me ~" H °\^). (5.3) 

n— >-±oo 

The advantage of this representation is that because H > 0, the spectrum of (e _/3ff ) is in the interval [0, 1]. For large 
fixed n, e~ me can be uniformly approximated by a polynomial in e~@ H . The advantage is that powers of e~@ H 
can be computed directly using the generating functional 

N b N c 

(C\e^ nH \B) =J^Y1 b ^Z[ 9l ^ mik - 9/,] (5.4) 

3 = 1 k=l 

without using analytic continuation. 

This suggest the following sequence of approximations to compute scattering amplitudes. First use narrow wave 
packets sharply peaked in linear momentum to approximate sharp momentum transition matrix elements in terms of 
S matrix elements in normalizablc states: 

( Pl! Ml , • • • , Pw , MjTlPi, A*, Pa, A*> « a^iJCE^rp^ ' ( 5 

Next approximate = (\Pj + |\E , /_) using (|5 . 3[) for large enough n. This step also involves solving the one-body 

problem for each asymptotic particle in the initial and final states: 

Next, after fixing n, approximate e l2nx on x G [0, 1] by a polynomial in x, which gives 

e Mne-^ ^J2cm(n)(e-^ H ). (5.7) 

Taken together these approximations, when performed in the correct order, provide a means to compute on-shell 
transition matrix elements using purely Euclidean methods. 



VI. TEST OF APPROXIMATIONS 

A mathematically controlled approximation is not automatically useful in all applications. To test the suggested 
sequence of approximations at the relevant GeV scale we consider a simple model based on a separable potential 



This is an exactly solvable model; a first test of the proposed method is to calculate scattering amplitudes in this 
model using matrix elements of e~^ H in normalizable states. 

Table 1: Degree 300 polynomial compared to e~ mx , n = 220 
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RealS-l vs n-limit - 1 GeV 
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FIG. 1: Real part of S compared to Kato 
Birman approximation as function of n. 



Real part of <klT(k)lk> (exact - black, polynomial - red) 




FIG. 3: Calculated compared to approximate 
calculations for the real part of transition 
matrix element for different energies. 



Imaginary S vs n-limit - 1 GeV 
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FIG. 2: Imaginary part of S compared to 
Kato Birman approximation as function of 
n. 

Im part of <klT(k)lk> (exact - black, polynomial - red) 




FIG. 4: Calculated compared to approximate 
calculations for the real part of transition 
matrix element for different energies. 



In this model sharp T matrix elements can be calculated with an error of approximately 1% using wave packets 
whose momentum widths are about 1/10 of the cm momentum. Figures 1 and 2 show the convergence of the real 
and imaginary parts of the S matrix evaluated in these wave packets as a function of n in equation (|5.6[) . Values of 
n between 200-300 are adequate in this model. /3 is a parameter that can be adjusted to improve the convergence. 
Polynomial approximations to e 2me are performed using Chcbyshcv expansions: 

N N 

e*» « icbToOr) + $>T fe (z) c, = ^- £ e **««* » costf^^). (6.1) 

k=l fc=l 

Polynomials of degree approximately 300 agree with e 2me P uniformly to better than 13 significant figures. 

Typical results are shown in table 1. Figures 3 and 4 compare the exact value of the real and imaginary parts of 
the sharp momentum transition matrix to the calculated values for momenta up to 2 GeV. For most values of k the 
exact and approximate values cannot be distinguished. 

These results suggest that it may be feasible to use this method to formulate relativistic few-body models. The open 
problems that have not been addressed in this preliminary work involve finding model Green functions or generating 
functionals satisfying the required conditions. Reflection positivity appears to be a fairly restrictive condition that 
requires additional study. The toy model discussed above did not require solutions of the one-body problem. How 
approximate solutions of the one-body problem are affected by the other approximations also requires further study. 

This work supported in part by the U.S. Department of Energy, under contract DE-FG02-86ER40286. 
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